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FRACTURE  MECHANICS  OF  FUNCTIONALLY 
GRADED  MATERIALS 


ABSTRACT 

Functionally  graded  materials  are  generally  two-phase  composites  with 
continuously  varying  volume  fractions.  Used  as  coatings  and  interfacial  zones,  they  help 
to  reduce  mechanically  and  thermally  induced  stresses  caused  by  the  material  property 
mismatch  and  to  improve  the  bonding  strength.  In  this  project  some  basic  problems 
concerning  fracture  mechanics  of  graded  materials  are  identified,  general  analytical 
methods  for  solving  the  related  crack  problems  are  developed,  the  singular  behavior  of 
the  solutions  for  typical  material  nonhomogeneities  is  examined,  and  solutions  of  some 
benchmark  problems  are  obtained.  The  results  are  intended  to  provide  technical  support 
for  material  scientists  and  engineers  who  are  trying  to  develop  methods  for  processing 
these  materials  and  for  design  engineers  who  are  interested  in  using  them  in  technological 
applications.  Typical  applications  of  fimctionally  graded  materials  include  thermal 
barrier  coatings  of  high  temperature  components  in  gas  turbines,  surface  hardening  for 
tribological  protection,  and  as  interlayers  in  microelectronic  and  optoelectronic 
components. 

The  results  foimd  show  that  by  eliminating  the  discontinuities  in  material  property 
distributions  the  mathemaitcal  anomalies  regarding  the  crack  tip  stress  oscillations  for  the 
interface  cracks  and  the  nonsquare-root  singularities  for  cracks  intersecting  the  interfaces 
are  also  eliminated.  From  the  viewpoint  of  fracture  mechanics  the  importance  of  this 
result  lies  in  the  fact  that  in  analyzing  the  components  involving  fimctionally  graded 
materials  one  can  use  the  existing  crack  tip  finite  element  modeling  developed  for 
ordinary  square-root  singularities  and  apply  the  energy  balance-based  theories  of 
conventional  fracture  mechanics.  The  benchmark  problems  considered  in  this  project  are 
concerned  primarily  with  surface  cracking  and  debonding. 
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1.  INTRODUCTION  TO  FUNCTIONALLY  GRADED  MATERIALS 

To  meet  the  expected  stringent  demands  of  future  technologies  in  power 
generation,  transportation,  aerospace,  and  microelectronics,  in  current  research  a  greater 
emphasis  will  have  to  be  placed  on  material  design,  more  specifically,  on  developing  new 
materials  or  material  systems  tailored  for  specific  applications.  Generally  such  materials 
tend  to  be  composites  and  intermetallics  having  homogeneous  bulk  properties.  The 
composites  may  be  fiber  or  filament  reinforced,  particulate,  or  layered  in  structure.  Many 
of  the  laminated  materials,  thin  films  and  coatings  fall  into  the  latter  category.  A 
common  feature  of  composites  is  that  they  consist  of  bonded  dissimilar  homogeneous 
materials.  Consequently,  in  studying  mechanics,  particularly  failure  mechanics  of  such 
materials  the  nature  of  interfaces  or  interfacial  regions  would  play  an  extremely  important 
role.  From  a  mechanics  view  point,  material  property  discontinuities  across  the  interfaces 
would  have  generally  two  undesirable  consequences,  namely  higher  residual,  thermal  and 
mechanical  stresses,  and  weaker  bonding  strength.  To  circumvent  these  difficulties,  very 
often  the  interfacial  regions  are  modified  by  intorducing  a  third  medium  in  the  form  of  an 
interlayer,  mechanically  roughening  the  surfaces,  or  using  coupling  agents. 

An  alternative  concept  which  may  also  be  used  to  overcome  some  of  the 
shorteomings  of  piecewise  homogeneous  materials,  partieularly  of  layered  dissimilar 
materials,  would  be  the  introduction  of  interfacial  regions  or  coatings  with  graded 
thermomechanical  properties'*^  [l]-[7].  Thus,  by  varying  the  volume  fractions  of  the 
constituents  between  zero  to  one  hundred  percent,  thereby  obtaining  a  continuous 
through-thickness  material  property  variation,  it  is  possible  to  obtain  not  only  smoother 
stress  distributions  and  lower  stress  concentrations  [8*],  [9],  but  also  higher  bonding 
strength  [10].  For  example,  in  [8*]  it  was  shown  that  the  point  of  intersection  of  the  free 
end  and  the  interface  between  a  homogeneous  coating  and  the  substrate  is  a  point  of 
stress  singularity  and,  consequently,  a  potential  location  of  debonding  crack  initiation. 


^*^  Numbers  in  brackets  refer  to  the  references  listed  in  Section  4  of  this  report.  The 
references  marked  by  an  (*)  describe  research  which  is  fully  or  partly  supported  by  the 
current  AFOSR  grant. 
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On  the  other  hand,  if  the  homogeneous  coating  is  replaced  by  a  graded  layer,  the 
singularity  disappears  and  the  stress  distribution  becomes  considerably  smoother. 

In  many  of  the  deposition  and  bonding  processes  used  in  ceramic  coatings,  it  is 
difficult  to  obtain  the  desired  strength  for  the  interfacial  region.  This  is  due  largely  to 
poor  adhesion  and  partly  to  high  stress  concentrations.  The  adverse  influence  of  both  of 
these  factors  can  be  reduced  significantly  by  introducing  a  graded  interfacial  zone 
between  the  two  materials.  The  technique  can  be  particularly  useful  for  material  pairs  in 
which  the  bonding  is  inherently  difficult,  for  example,  in  [10]  it  was  shown  that  for  a 
diamond  film  synthesized  over  a  50/50  W/Mo  alloy  by  using  a  DC  plasma  jet,  the 
measured  bonding  strength  is  less  than  10  kg/cm  .  On  the  other  hand,  if  a  graded 
interfacial  zone  is  introduced  by  first  plasma-spraying  the  substrate  with  tungsten  carbide 
and  then  gradually  adding  increasing  amounts  of  methane  and  hydrogen  before  growing 
the  diamond  film,  the  adhesive  strength  is  measured  to  be  over  150  kg/cm  . 

These  particulate  composites  with  continuously  varying  volume  fractions  are 
called  fiinctionally  graded  materials  (FGM).  By  controlling  not  only  the  composition  of 
profile  but  also  the  microstructure,  the  concept  of  FGM’s  could  provide  a  great  deal  of 
flexibility  in  material  design.  As  the  processing  techniques  improve,  the  potential  for 
special  applications  of  FGM’s  appears  to  be  nearly  unlimited.  However,  in  the 
immediate  future  the  primary  use  of  these  new  materials  will  most  likely  be  limited  to 
thermal  barrier  coatings,  wear-resistant  coatings  and  corrosion-resistant  coatings  [11]. 

1.1  Thermal  Barrier  Coatings 

Within  the  past  decade  or  so  considerable  progress  has  been  made  in  using 
ceramic  coatings  to  protect  metallic  components  from  high  temperature  turbine 
environments.  These  thermal  barrier  coatings  (TBC’s)  are  currently  being  used  in 
conjunction  with  air  cooling  to  prolong  the  life  of  hot-section  turbine  components  in 
aircraft  engines.  The  application  of  TBC’s  also  offers  the  possibility  of  increasing  the 
thermodynamic  efficiency  of  land-based  turbines  by  increasing  the  inlet  temperature  of 
gases.  Thus,  the  TBC  technology  is  considered  by  many  as  a  viable  means  for 
developing  more  efficient  aircraft  engines  and  stationary  gas  turbines.  There  are. 
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however,  several  major  technical  issues  involving  the  next  generation  of  TBC’s,  namely 
(i)  improvements  in  processing  techniques  from  both  economic  and  performance 
standpoints,  (ii)  understanding  the  failure  mechanisms  of  TBC’s  in  simulated  and  actual 
turbine  environments,  and  (iii)  developing  characterization  and  test  methodologies  to 
measure  time  and  temperature-dependent  characteristics  of  TBC’s  for  reliable  life 
prediction. 

The  current  approach  for  accommodating  the  material  property  mismatch  between 
the  ceramic  coating  and  metallic  substrate  is  to  make  the  ceramic  layer  to  be  more 
compliant  (or  “strain  tolerant”)  by  introducing  defects  such  as  microcracks,  porosity  and 
segmented  columnar  microstructure.  The  difficulty  with  this  solution  is  that  the 
particular  microstructural  features  that  provide  compliant  coatings  also  provide  rapid 
diffusion  paths  for  oxygen.  The  experience  seems  to  indicate  that  by  far  the  most  critical 
factor  limiting  the  performance  of  the  state-of-the-art  TBC’s  is  the  spallation  of  the 
ceramic  layer  which  takes  place  either  along  a  plane  parallel  to  the  ceramic-bond  coat 
interface  (as  in  plasma-spray  coatings),  or  along  the  interface  between  the  bond  coat  and 
its  oxide  scale  (usually  A^203)  (as  in  coatings  processed  by  using  electron  beam  physical 
vapor  deposition).  In  both  cases  the  life  of  the  coating  is  controlled  by  the  oxidation  time 
and  the  number  of  thermal  cycles. 

The  basic  conjecture  behind  the  FGM  concept  is  that  by  replacing  sharp  interfaces 
with  graded  interfacial  zones  or  by  replacing  homogeneous  ceramic  layers  with 
metal/ceramic  FGM’s,  it  is  possible  to  improve  the  resistance  of  the  coating  to  spallation 
as  a  result  of  reduced  stress  levels  and  improved  bonding  strength. 

1.2  Wear-Resistant  Coatings/Abradable  Seals 

One  obvious  application  of  ceramic  coatings  apears  to  be  to  provide  the  necessary 
hardness  or  wear  resistance  to  the  surfaces  of  structural  components  transmitting  forces 
through  contact,  such  as  gears,  bearings,  cams  and  machine  tools.  Intuitively  it  is  clear 
that  the  fatigue  life  of  these  components  may  be  improved  quite  considerably  if  one  uses 
a  graded  rather  than  a  homogeneous  ceramic  coating  on  the  metallic  substrate.  In  these 
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load-transfer  components  FGM  coatings  would  provide  the  necessary  surface  hardness 
without  sacrificing  the  overall  toughness. 

A  related  application  of  the  FGM  interfacial  zones  or  coatings  may  be  found  in 
the  abradable  seals  used  in  some  stationary  gas  turbines  to  help  minimize  the  gas  leakage 
through  the  gap  between  the  tips  of  rotating  blades  and  the  turbine  shroud.  The  problem 
is  again  a  contact  problem  and  the  primary  material  property  requirements  are 
abradability  and  toughness.  The  main  components  in  the  shroud  are  the  metallic  structure 
or  the  substrate,  bond  coat,  high  density  ceramic  and  very  low  density  ceramic  with  a 
graded  zone  replacing  every  sharp  interface  [11]. 

In  the  past  wear  and  corrosion-resistant  coatings  have  been  used  quite  extensively 
in  industrial  machinery.  The  coating  material  has  been  metals  such  as  stainless  steels. 
Mo-based  alloys  and  WC-Co  as  well  as  ceramics  such  as  A?203-Ti02  and  Cr203.  For 
example,  WC-Co  and  Cr3C2  /  NiCr  have  been  extensively  used  in  aircraft  industry  to  coat 
various  turbine/compressor  components  and  mid-span  stiffeners  for  improved  wear- 
resistance.  Other  applications  of  wear-resistant  coatings  have  been  in  printing  rolls,  steel 
mills,  petrochemical  industry,  and  automobile  industry.  Most  of  these  coatings  have  been 
deposited  by  using  a  thermal  spray  process.  Since  thermal  spray  technique  is  readily 
suitable  for  grading  the  composition  of  the  coating,  service  life  improvements  can  be 
obtained  in  all  applications  of  wear-resistant  coatings  by  using  the  FGM  concept. 

1.3  Oxidation-Resistant  Coatings 

The  main  desirable  characteristics  of  an  ideal  TBC  appear  to  be  low  conductivity 
for  thermal  insulation,  high  coefficient  of  thermal  expansion  to  match  that  of  the  metallic 
substrate,  and  high  resistance  to  oxygen  diffusion.  It  is  highly  unlikely  to  find  all  these 
favorable  properties  in  a  single  material.  Many  of  the  well  known  ceramics  have  either 
high  conductivity  and  low  oxygen  diffusivity  or  low  conductivity  and  high  oxygen 
diffiisivity.  to  prevent  oxygen  diffusion  at  some  point  a  layer  of  AI2O3  or  mullite 
(A203»2Si02)  may  be  needed.  However,  these  materials  have  considerably  higher 
thermal  conductivity  than  that  of,  for  example,  YSZ  (ytria-stabilized  zirconia).  Thus,  the 
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problem  appears  to  be  an  optimal  design  of  a  multi-layered  structure,  including  graded 
interfacial  zones  and  coatings. 

Typically,  the  current  design  of  TBC’s  consists  of  a  partially  stabilized  zirconia 
coating  deposited  on  an  intermediate  metallic  bond  coat  (e.g.,  NiCr  AlY)  which  is  plasma 
sprayed  on  the  (superalloy)  substrate  [7].  the  main  function  of  the  bond  coat  is  to  protect 
the  substrate  against  oxidation.  It  also  helps  to  reduce  thermal  expansion  mismatch 
between  the  ceramic  coating  and  the  metallic  substrate,  and  provides  surface  texture  to 
improve  bonding.  At  high  temperatures  an  oxide  (AJ203)  scale  is  formed  along  the  PSZ- 
bond  coat  interface.  Even  though  this  AI2O3  layer  forms  an  oxygen  diffusion  barrier,  it 
also  introduces  a  weak  cleavage  plane  which,  under  thermal  cycling,  may  lead  to 
spallation.  This  difficulty  may  be  overcome  by  introducing  a  graded  (NiCrAf2Y-PSZ) 
layer  between  the  bond  coat  and  the  ceramic  layer  [7]. 

2.  OBJECTIVES  OF  THE  RESEARCH  PROGRAM 

The  objectives  of  the  research  program  supported  by  the  AFOSR  Grant  have  been 

•  To  identify  some  basic  problems  concerning  the  fracture  mechanics  of 
functionally  graded  materials. 

•  To  develop  general  methods  for  solving  the  related  crack  problems. 

•  To  examine  the  singular  behavior  of  the  solutions  of  the  crack  problems  in 
FGM’s  and  compare  them  with  the  corresponding  results  in  homogeneous  and 
piecewise  homogeneous  materials. 

•  To  obtain  solutions  for  some  benchmark  problems. 

In  general,  the  results  of  the  research  program  are  intended  to  provide  technical 
support  for  material  scientists  and  engineers  who  are  trying  to  develop  methods  for 
processing  FGMs  and  for  design  engineers  who  are  interested  in  applications. 
Specifically,  the  crack  tip  singularities  foxmd  and  the  benchmark  solutions  provided 
would  be  quite  useful  in  the  development  and  testing  of  finite  element  moldes  for  solving 
more  complex  problems  involving  FGMs. 
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3.  FRACTURE  MECHANICS  OF  GRADED  MATERIALS 

-  STRESS  SINGULARITIES  AND  BENCHMARK  SOLUTIONS 

In  this  section,  after  describing  some  elementary  concepts  of  fracture  mechanics, 
a  brief  discussion  of  the  topic  as  it  relates  to  the  graded  materials  is  presented.  Stress 
singularities  for  each  typical  crack  geometry  are  then  examined  and  sample  results 
obtained  from  a  series  of  benchmark  solutions  are  discussed,  the  analytical  details  and 
more  extensive  results  are  given  in  References  [8  ]  and  [12*]-[20  ]. 

3.1  Fracture  Mechanics  -  Elementary  Concepts 

In  a  broad  sense  “fracture”  is  creation  of  new  surfaces  in  solids.  The  macroscopic 
theories  of  fracture  are  based  on  the  notions  of  continuum  mechanics  and  classical 
thermodynamics.  In  studying  the  subject,  it  is  implicitly  assumed  that  the  material 
contains  some  macroscopic  flaws  and  the  medium  is  homogeneous  in  the  sense  that  the 
size  of  these  flaws  is  large  compared  to  the  characteristic  microstructural  dimension  of 
the  materials.  The  problem  is  to  study  the  effect  of  the  applied  loads,  flaw  geometry  and 
the  environmental  conditions  on  the  fracture  of  the  solid  -  a  subject  which  has  come  to  be 
known  as  fracture  mechanics.  The  fundamental  criterion  of  fracture  initiation  and 
propagation  is  based  on  the  energy  balance  concept.  Let  the  solid  contain  a  dominant 
flaw  which  is  usually  considered  to  be  a  planar  crack  of  surface  area  A.  Under  given 
external  loads  if  the  crack  grows  by  an  amoxmt  dA  in  time  dt,  the  thermodynamic 
equilibrium  of  the  solid  requires  that 

dU_dV  dT  dP 
dtdtdtdt’ 

where  U,  V,  T  and  D  respectively  are  the  work  of  the  external  loads,  the  recoverable 
internal  energy,  the  kinetic  energy,  and  the  sum  of  all  dissipated  energies  such  as  surface 
tension,  plastic  work,  viscous  dissipation,  etc.  If  the  energy  dissipation  takes  place  only 
around  the  advancing  periphery  of  the  crack,  in  the  quasi-static  case  T  is  negligible  and 
defining  dD/dA=Ac>  (1)  may  be  expressed  as 

(2) 
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In  the  fracture  criterion  given  by  (2)  the  left  hand  side  is  the  energy  available  and 

Ac  is  the  energy  required  to  create  a  unit  area  of  new  fracture  surface.  They  are  also 

known  as  the  crack  driving  force  and  the  fracture  toughness,  respectively. 

If  the  solid  is  brittle  or  quasi  brittle,  the  characteristic  size  of  the  zone  of  energy 
dissipation  around  the  crack  front  is  generally  small  in  comparison  with  the  length 
parameter  of  the  crack.  It  may,  therefore,  be  assumed  that  the  energy  (U-V)  flowing  into 
the  fractme  zone  will  come  from  the  elastic  bulk  of  the  solid  and  will  not  be  critically 
dependent  on  the  details  of  the  stress  state  very  near  the  crack  tip.  Furthermore,  it  is  also 
clear  that  in  this  case  the  stress  state  in  the  elastic  bulk  of  the  actual  solid  will  not  differ 
from  a  purely  elastic  solution  to  any  significant  extent.  Consequently,  one  may  be 
justified  in  calculating  the  crack  driving  force  by  using  the  linear  elastic  crack  solution. 
By  using  the  concept  of  crack  closure  it  can  then  be  shown  that  the  increment  d(U-V)  of 
the  energy  available  for  fracture  may  be  evaluated  from  the  asymptotic  stresses  and  the 
crack  opening  displacements  near  the  crack  tip  which,  in  turn,  may  be  obtained  from  the 
three  dimensional  elasticity  solution  as  follows: 


o,(r,0)s^/,„(e)+^/,,(e) ,  (;,y  =  x,>.)  , 
Os('',e)s-^/3,(6)  ,  0  =  x,y)  , 

E  E 

w"'  -  w"  s  — 


(3) 

(4) 

(5) 

(6) 


where  and  k^  are  the  modes  I,  II,  and  III  stress  intensity  factors,  fy ,  f-^y  and  /j, 

are  known  functions  and  E,  v  and  p  are  the  elastic  constants,  E  =  2p(l  +  v) .  From  the 
crack  closure  energy  it  may  then  be  shown  that 


-k{  ,  *2  = 


E 


U2  ^ 

2p 


oA 


(7) 

(8) 
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where  b  is  the  total  energy  available  for  fracture. 

Equation  (7)  indicates  that  one  may  also  use  k,  in  place  of  b^  as  the  measure  of 
the  crack  driving  force.  For  mode  I  loading  conditions,  for  example,  defining 

Ki=k^^  ,  A,,  =  bjc  ,  K,c  =  »  (9) 

the  fracture  criterion  (2)  may  be  expressed  as 

K,<K,c  .  (10) 

Equation  (10)  has  proved  to  be  very  useful  in  considering  the  fracture  stability. 
However,  perhaps  the  most  useful  application  of  the  stress  intensity  factors  may  be  found 
in  analyzing  the  subcritical  crack  growth  processes. 

3.2  Fracture  Mechanics  of  Graded  Materials-Embedded  Cracks 

In  studying  the  fracture  mechanics  of  FGMs  one  may  have  to  deal  with  a  number 
of  distinct  problem  areas.  The  first  is  the  investigation  of  the  nature  of  stress  singularities 
near  the  tip  of  a  crack  embedded  ina  nonhomogeneous  medium.  The  second  is  the 
general  problem  of  debonding  and  the  effect  of  a  possible  “kink”  in  material  property 
distributions  on  the  behavior  of  stress  singularities.  The  third  is  the  basic  surface 
cracking  problem  and  the  nature  of  the  stress  singularities  for  cracks  intersecting  the 
interfaces.  The  fourth  is  the  problem  of  stress  relaxation  due  to  periodic  surface  cracking 
in  FGM  coatings.  The  fifth  significant  problem  area  is  concerned  with  the  stress 
concentrations  and  the  resulting  delamination  cracking  at  the  free  ends  of  FGM  coated 
homogeneous  substrates  imder  residual  or  thermal  stresses.  In  the  following  sections 
these  problem  areas  are  briefly  discussed  and  some  examples  are  given. 

To  examine  the  influence  of  the  material  nonhomogeneity  on  the  asymptotic 
stress  state  near  the  crack  tips,  we  first  consider  the  plane  elasticity  problem  for  an 
infinite  medium  containing  a  line  crack  (Fig.  1).  For  simplicity  we  will  assume  that  the 
Poisson’s  ratio  v  of  the  medium  is  constant  and  the  shear  modulus  is  approximated  by 
p(x,y)  =  Poexp(px+yv)  ,  (11) 
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where  iXqj  P  andy  are  known  constants.  This  problem  was  solved  for  a  crack  along 
y=0,  -a<x<a  under  arbitrary  loading  conditions  [12*],  [22*].  It  was  shown  that  near  the 
crack  tip  x=a  the  stresses  have  the  following  asymptotic  behavior: 


a  y(x,y)  =  exp[r(P  cos0  +y  sin0)] 


,  (i,j  =  x,y), 


where  the  stress  intensity  factors  it,  and  k2  are  defined  by 

k.  (a)  =  limV2(x - a)a  ^(x,0)  ,  (a)  =  lim V2(x  - a)a  (x,0)  , 


(12) 

(13) 


and  the  functions  f^y  and  are  identical  to  those  found  for  the  homogeneous  materials 

given  in  (3).  Note  that  the  asymptotic  stress  states  for  homogeneous  materials  (3)  and 
FGMs  (12)  are  identical  only  at  r=0.  However,  since  the  crack  opening  displacement  is 
also  influenced  in  a  way  similar  to  stresses,  the  crack  driving  force  (or,  for  “fixed  grip” 
conditions,  the  strain  energy  release  rate)  was  found  to  be  identical  to  that  calculated  for 
the  homogeneous  materials,  namely 


8p(a,0) 


(14) 


where  k=3-4v  for  plane  strain  and  k=(3-v)/(1+v)  for  plane  stress  conditions. 

Some  sample  results  for  an  embedded  crack  of  length  2a  in  an  infinite  medium 
imder  plane  strain  conditions  are  given  in  Table  1.  Referring  to  Fig.  1  and  (1 1)  we  first 
define 


x'  =  x  cos  0  +  y  sin  0  ,  p  =  to  cos  0 ,  y  =  to  sin  0 ,  to^  =  p^  +y  (15) 

=  M(^')  =  lio  exp(tox')  .  (16) 

The  medium  is  assumed  to  be  loaded  by  fixed  grips  away  from  the  crack  region  with 
8^,y(jc',+oo)  =  Eq.  Thus,  the  normalizing  stress  intensity  factor  for  the  results  given  in 
the  table  is 

^0  =2(l+v)poeoVa  -  (17) 

Note  that  the  crack  orientation  angle  0=0  corresponds  to  a  mode  I  problem  for  which 
(+a)  =  0 ,  whereas  for  0  =  7t  /  2  the  loading  is  parallel  to  the  crack  and,  consequently. 
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Table  1 .  The  effect  of  material  nonhomogeneity  parameter  a©  and  the  crack  orientation 
angle  0  on  the  stress  intensity  factors;  v=0.3;  loading:  uniform  strain  away  from  the  crack 
region,  s  yy  (x',+oo)  =  s  „ . 

aco 

0/7t 

ki(a)/ko 

ki(-a)/ko 

k2(a)/ko 

k2(-a)/ko 

0 

1.196 

0.825 

0 

0 

0.1 

1.081 

0.750 

-0.321 

-0.254 

0.2 

0.781 

0.548 

-0.514 

-0.422 

0.25 

0.3 

0.414 

0.290 

-0.504 

-0.437 

0.4 

0.121 

0.075 

-0.304 

-0.282 

0.5 

0 

0 

0 

0 

0 

1.424 

0.674 

0 

0 

0.1 

1.285 

0.617 

-0.344 

-0.213 

0.2 

0.925 

0.460 

-0.548 

-0.365 

0.5 

0.3 

0.490 

0.247 

-0.532 

-0.397 

0.4 

0.146 

0.059 

-0.314 

-0.269 

0.5 

0 

0 

0 

0 

0 

6.317 

0 

0 

0.1 

5.376 

-0.037 

0.2 

3.315 

0.115 

-1.155 

-0.090 

2.5 

0.3 

1.441 

0.082 

-0.900 

-0.158 

0.4 

0.369 

0.004 

-0.429 

-0.179 

0.5 

0 

0 

0 

0 

Table  2.  Stress  intensity  factors  for  a  plane  strain  and  a  permy-shaped  crack 
medium  under  tension  Gq  perpendicular  to  the  plane  of  the  crack;  v=0.3. 

in  a  graded 

©a  0 

0.1 

0.25 

0.5 

1.0 

2.5 

5.0 

Plane  Strain  Crack 

kj/CToV^  1 

1.008 

0.036 

1.101 

1.258 

1.808 

2.869 

k2/C7o7o  0 

0.026 

0.065 

0.129 

0.263 

0.697 

1.567 

Penny-Shaped  Crack 

1.002 

1.012 

1.038 

1.118 

1.442 

2.083 

0.017 

0.041 

0.083 

0.168 

0.440 

0.960 

all  stress  intensity  factors  are  zero. 


11 


Table  2  shows  some  limited  results  comparing  plane  strain  and  penny-shaped 
crack  solutions  for  a  graded  material  under  uniform  tension  (Tq  perpendicular  to  the  plane 
of  the  crack  where  co=0  corresponds  to  the  results  for  a  homogeneous  medium.  The  table 
shows  that  the  stress  intensity  factors  in  FGMs  are  higher  than  that  in  homogeneous 
materials  and  that  the  influence  of  the  material  nonhomogeneity  on  the  stress  intensity 
factors  is  more  severe  for  a  plane  strain  than  for  a  penny-shaped  crack.  For  further  results 
see  [12],  [22],  [23]  and  [24].  Also,  for  the  solution  of  embedded  crack  problems  in  an 
FGM  layer  under  mechanical  or  thermal  loading  see  [15*]  and  [16*]. 


3.3  The  Debonding  Problems  for  FGMs 

Consider  the  crack  problems  shown  in  Figs.  2a  and  2b.  Figure  2a  describes  the 
debonding  problem  in  piecewise  homogeneous  materials,  whereas  Fig.  2b  refers  to  a 
FGM  bonded  to  a  homogeneous  substrate.  In  both  cases  h=0  refers  to  an  “interface 
crack”.  In  terms  of  the  unknown  functions 


/iW  =  ;^(v^-V)  , /2(x)  =  ^K-m’)  , 
dx  dx 


(18) 


in  each  case  the  formulation  of  the  problem  may  be  reduced  to  a  system  of  integral 
equations  of  the  form 

1 

1 


71  •!.  ,  t-x  2|ii(0) 


(/■  =  1,2)  ,  -a  <x  <a  (19) 

where  the  kernels  are  known  functions  which  depend  on  h  and  the  material 
parameters,  k-j  is  associated  with  the  infinite  medium,  k{j  represents  the  geometry  of  the 
medium,  and 

P\  i^)=^yy  ’P2(x)=<yxy  (20) 

are  the  crack  surface  tractions  which  may  be  expressed  in  terms  of  the  external  load.  The 
kernels  k(j  are  bounded  for  all  values  of  h.  For  h>0  the  functions  k-j  are  eiIso  bounded. 

Thus,  for  h>0  the  crack  is  an  embedded  crack  and  (19)  would  lead  to  the  asymptotic 
stresses  given  by  (3)  and  (12)  for  problems  described  by  Figs.  2a  and  2b,  respectively. 
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On  the  other  hand,  for  h=0  in  problem  2a  the  kernels  k^^  and  Atjj  would  become  a 

Cauchy  kernel  (t-xy^  and  ^21  would  degenerate  to  a  delta  function  8(t-jc) 

[25],  [26].  Consequently,  in  this  interface  crack  problem  the  integral  equations  become 
one  of  the  second  kind  leading  to  the  well-known  anomalous  stress  oscillation  behavior 
very  near  the  crack  tips. 


For  h=0  in  problem  2b,  however,  the  leading  terms  of  the  kernels  become 


=  - 

— 1^22 


ny 


\t-x 


8  t-x 


Ki  =  -^21  =  +  jlogk  “  ^1  ’  Y  =  t3n(j) 


0  » 


(21) 


which  would  indicate  that  (19)  would  remain  to  be  an  ordinary  system  of  singular  integral 
equations  of  the  first  kind  and  would  have  the  asymptotic  solution  given  by  (12).  It  is, 
therefore,  seen  that  the  anomalous  behavior  of  the  crack  tip  stress  oscillations  may  be 
eliminated  by  “smoothing”  the  material  property  distribution  (or  by  removing  the 
property  discontinuity).  A  qualitative  description  of  the  interface  crack  geometries  and 


the  singular  kernels  ky  may  be  seen  in  Fig.  3. 

Figure  4  describes  an  interface  crack  problem  for  a  FGM  coating  bonded  to  a 
homogeneous  substrate.  In  actual  applications  the  medium  contains  an  additional  layer  of 
bond  coat  between  the  substrate  and  the  FGM  coating.  Typically  the  materials  for  the 
bond  coat  and  the  substrate  are  NiCrAfY  and  MCrA^Y,  respectively,  wehre  M  stands  for  a 
“metal”  (e.g.,  NiCo).  Since  in  most  cases  these  two  materials  have  approximately  the 
same  thermomechanical  constants,  in  the  example  considered  the  bond  coat  is  assumed  to 
be  an  extension  of  the  substrate.  Because  of  the  oxide  scales  (usually  AI203)  forming 
along  bond  coat/FGM  interface,  it  is  also  assumed  that  this  interface  is  the  weak  cleavage 
plane  and  contains  a  crack  of  length  2a.  In  the  example  h^  and  h2  are  the  thicknesses  and 
Pi  and  P2(y)  are  the  shear  moduli  of  the  substrate  and  the  FGM,  respectively.  The 
nonhomogeneity  constant  y  is  defined  by 

li2 (t)  =  111  exp(yv)  ,  y  =  ■^log[p2 (/i^)  /  p, ]  .  (22) 
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The  problem  is  solved  under  uniform  crack  surface  pressure,  a  ^  (x,0)  =  -<j  ^  or  shear, 


G  ^y(x,0)  =  -Tq  .  In  each  case  the  calculated  stress  intensity  factors  kj,  k2  and  the  strain 
energy  release  rate  t  are  normalized  with  respect  to  the  following; 

Kq  g  q  V n  ,  Kq  —  ^  0  V n  ,  Aq  —  cj  gUt ,  Ag  —  X  qQ  ,  (23) 

8p,  8p, 

which  correspond  to  the  values  obtained  for  an  infinite  medium  under  plane  strain 
conditions  and  imiform  tractions  a  ^  (x,+oo)  =  a  g  or  <t  =  x  g . 

Some  sample  results  showing  the  effect  of  material  nonhomogeneity  constant  y 
and  the  relative  dimensions  hj/a  and  h2/a  are  given  in  Figures  5-10.  In  all  cases  y=0 
corresponds  to  the  homogeneous  medium.  As  physically  expected,  as  y  increases,  the 
stiffness  of  the  coating  increases,  leading  to  decreasing  crack  driving  force  h. 

Clearly  the  problem  is  one  of  mixed-mode  and  unless  the  plane  of  the  crack  is  a 
weak  cleavage  plane  (e.g.,  due  to  oxidation),  the  crack  may  not  grow  in  a  co-planar 
fashion.  In  this  case,  the  stress  intensity  factors  kj  and  k2  may  be  used  to  determine  the 
probable  crack  growth  angle,  Oq.  This  angle,  measured  from  the  positive  x  axis  may  be 
calculated  from  gg  /  59  =0,  a  gg  (r,0  g)  >  0  [27],  or 

A:,  sinBg  +  A:2(3cos0g  -1)  =  0  , 


A:,  cos^  ^ ~ sin0g  >  0  .  (24) 

Some  calculated  values  of  0o  for  v=0.3,  hi/a=100,  h2/a=0.25, 0.5, 1, 2, 10  and  for  uniform 
crack  surface  presure  CTq  are  shown  in  Fig.  1 1 .  The  figure  indicates  that  if  the  medium  is 
isotropic  with  regard  to  the  crack  growth  resistance  Ag  near  the  crack  tip,  maximum  strain 
energy  release  and,  as  a  result,  further  crack  growth  would  take  place  in  a  direction 
toward  the  less  stiff  material.  On  the  other  hand,  if  is  0-dependent,  then  the  crack 
growth  direction  0o  would  be  determined  by  maximizing  b(Q)/bj(Q),  b(Q)  being  the  strain 
energy  release  rate  for  a  small  radial  crack  extension  in  the  direction  of  0.  Further  results 
may  be  foirnd  in  [14*]. 

It  should  be  noted  that  in  formulating  this  problem  the  crack  surface  tractions  are 
assumed  to  be  arbitrary  functions  of  x.  Therefore,  the  technique  described  in  [14*]  may 
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be  used  to  solve  any  quasi-static  interface  crack  problem  involving  FGM  coatings  after 
properly  separating  the  solution  for  the  uncracked  problem.  Also,  for  the  validity  of  the 
solution  it  is  essential  that  the  normal  crack  opening  v^-v'  be  positive.  Otherwise  the 
problem  must  be  treated  as  a  crack/contact  problem,  the  condition  ki(-a)  determining  the 
crack  closure  point  [16  ]. 


3.4  Cracking  perpendicular  to  interfaces  and  surfaces 

In  ceramic  and  ceramic/metal  FGM  components  generally  a  common  mode  of 
failure  is  siuface  cracking  which  could  penetrate  to  the  interface  and  cause  spallation. 
The  main  problem  here  is  assessing  the  influence  of  the  material  nonhomogeneity  on  the 
fracture  mechanics  parameters  (such  as  b  and  kj)  for  surface  cracks  and  cracks 
terminating  at  an  interface.  Figures  2c  and  2d  show  the  crack  geometry  for  the  latter 
problem  in  piecewise  homogeneous  and  in  nonhomogeneous  materials.  Because  of 
symmetry,  generally  these  are  all  mode  I  problems.  Thus,  if  we  define  the  unknown 
function  and  the  crack  surface  traction  by 


g(x)  =  —  [v(x,+0) - v(x,-0)]  ,  p(x)  =  a  (x,0)  ,  a  <x<b  , 
ox 

the  integral  equation  for  the  general  problem  may  be  expressed  as 


—  \[-^  +  k^{x,t)  +  kAx,t)'\g{t)dt  =  ,  a<x<b  , 

7C  '  t-x  2\\.2 


(25) 


(26) 


where,  again,  k^  is  associated  with  two  bonded  semi-infinite  media,  kf  represents  the 
geometry  of  the  composite  medium,  and  kf  is  always  boxmded.  for  an  embedded  crack, 
a>0  and  k^  is  also  bounded.  However,  for  a-0,  ks  could  be  singular.  In  fact,  for  a=0  in 

piecewise  homogeneous  materials  (Fig.  2c)  it  is  known  that 

.2 


kXx,t)  = 


CjX 

'  •  +  — - — r-i- 


C3X 


,  (Q<{t,x)<b)  , 


(27) 


t+x  (t+xY  (t+xY 
where  Cj,  C2,  and  C3  are  bimaterial  constants  [28].  Note  that  as  t  and  x  approach  the  end 
point  x=0,  ks  tends  to  infinity  and,  hence,  would  contribute  to  the  singular  behavior  of  the 
solution  giving 
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where  gy  are  known  functions,  kj  is  a  “stress  intensity  factor”  and  the  power  of  stress 
singularity  a  is  a  fimction  of  the  bimaterial  constants.  Generally  a>l/2  for  and 
a<l/2  for  oc=1/2  being  the  value  for  P2'^Mi-  From  the  viewpoint  of  fracture 

mechanics,  the  consequence  of  having  is  that  as  the  crack  intersects  the  interface, 

the  stress  and  deformation  states  would  not  remain  self-similar  and,  hence,  it  would  not 
be  possible  to  use  the  fracture  theories  based  on  the  energy  balance  concept  to  calculate  a 
strain  energy  release  rate  or  to  use  the  stress  intensity  factors  as  the  crack  driving  force. 
This,  then,  is  the  second  anomalous  behavior  regarding  the  stress  state  near  the  crack  tip 
in  bonded  dissimilar  homogeneous  materials. 

If  we  now  “smooth”  the  material  property  distribution  and  assume  that  medium  1 
is  a  FGM  (Fig.  2d),  it  can  be  shown  that  for  a=0  the  leading  terms  of  k^  become  [29] 


k,{x,t)  = 


~  “I 


d.t 


d.tx 


■+d^  log(r  +  x) 


(29) 


t  +  x  t  +  x  {t  +  xf 

where  di...d4  are  bimaterial  constants.  Note  that  the  kernel  given  by  (29)  is  square 
integrable  and,  therefore,  would  have  no  contribution  to  the  stress  singularity  at  x=0. 
Consequently,  the  stresses  would  have  the  standard  square-root  singularity  and,  by 
smoothing  the  material  property  distribution  through  the  introduction  of  FGM,  the 
anomalous  behavior  of  the  stress  state  would  again  be  eliminated. 

Figure  12  shows  the  mode  I  stress  intensity  factor  for  a=0  and  p(x)=-CTo  in  Fig.  2d. 
The  normalized  stress  intensity  factors  shown  in  the  figure  are  defined  by 

k{^a)  =  k,{Qi)lcs^^IbT2  ,  k(b)  =  k,(b)lo^S72  ,  (30) 

k^ (0)  =  lim V-2x<t ,  (x,0)  ,  k^ib)  =  \im^j2(x~^a 2yyix,0)  .  (31) 

The  shear  modulus  of  FGM  in  Fig.  2d  is  assumed  to  be 

|Xi(x)  =  p2exp(px)  (32) 

where  p2  is  constant.  It  is  thus  seen  that  for  p^oo  pi^O  and  the  problem  becomes  an 


ordinary  edge  crack  problem  in  a  homogeneous  half  space  for  which 
A:,(0)^oo  ,  A:,(b)^li861<To>/h72  . 


(33) 


For  P=0  the  medium  is  homogeneous  and 

k,{Q))  =  k^(J})  =  o,-IbJl  .  (34) 

In  the  other  limiting  case  of  P=-oo,  becomes  infinite  and  for  the  resulting  problem  of  a 
crack  terminating  at  the  interface  we  have 

it, (0)^0,  A:, (6)^0.8710  CToV^  .  (35) 

The  analytical  details  and  further  results  for  this  problem  may  be  found  in  [29]. 

The  basic  problem  of  surface  cracking  in  a  FGM  layer  is  shown  in  Fig.  13.  In  this 
case  equations  (25)  defining  the  unknown  and  input  functions  and  (26)  describing  the 
form  of  the  integral  equation  are  still  valid.  In  the  surface  crack  problem  under 
consideration  it  can  be  shown  that  [15  ] 


k^{x,t)  = 


1 


6x 


4x^ 


(36) 


t  +  x  (t  +  xf  (t  +  xY 
which  is  identical  to  the  singular  kernel  found  for  a  homogeneous  medium  with  an  edge 
crack  [28]. 

Some  sample  results  for  the  layer  shown  in  Fig.  13  under  fixed-grip  loading  or 
constant  strain  e"  =8^,  membrane  loading  N  and  bending  moment  M  are  given  in 

Figures  14-19.  Figures  14-16  show  the  stress  component  Oyy  in  the  uncracked  layer.  In 
this  problem,  too,  it  is  assumed  that  the  Poisson’s  ratio  is  constant  and  the  Yoxmg’s 
modulus  at  x=0  and  x=h  is  given.by  E,  and  E2,  respectively,  E2/E1  being  the  measure  of 
material  nonhomogeneity.  The  Yoimg’s  modulus  E(x)  and  the  material  nonhomogeneity 
parameter  y  are  defined  by 


E(x)  =  ,  Y  =  jlog(E,  IE,).  (37) 

h 

The  normalizing  stresses  shown  in  the  figures  are  defined  by 

(To  =  £,So /(1-v^) ,  a,  =  ,  CTj  =  6M//i^  .  (38) 

Note  that  the  results  given  for  E2/Ei=l,  2,  5,  10,  20  in  Figures  15  and  16  may  also  be 
interpreted  for  E2/Ej=l,  0.5,  0.2,  0.1,  0.05,  respectively.  The  membrane  resultant  Nyy=N 
is  applied  along  x=h/2  away  from  the  crack  region.  Since  the  ends  are  free  to  rotate  the 
stress  distribution  shown  in  Fig.  15  exhibits  considerable  bending  effect. 
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The  mode  I  stress  intensity  factors  kj  corresponding  to  the  stresses  shown  in 
Figures  14-16  are  given  in  Figures  17-19,  respectively.  In  Figure  17,  the  normalizing 
stress  cTo  is  the  value  of  Cyy  at  x=0  in  the  uncracked  medium.  Thus,  for  all  values  of 

E2/E1,  as  the  crack  length  b  approaches  zero,  la^-Jb  would  approach  1.1215,  the  edge 
crack  solution  for  a  homogeneous  half  plane  [28].  Similarly,  note  that  under  membrane 
loading  and  bending,  as  b->0,  the  limiting  value  of  the  stress  intensity  factor  would  tend 
to 

=  1.1215  c^(0)^/6  ,  (6  «/2)  (38) 

where  ayy(x)  is  the  stress  in  the  rmcracked  medium.  It  should  be  observed  that  in  all 
cases  ki(b)  becomes  unbounded  as  b— >h  (Figures  17-19). 

Some  results  for  a  surface  crack  in  a  FGM  layer  under  thermal  loading  are  shown 
in  Figures  20-23.  Figure  20  shows  the  stress  distribution  CTyy(x)  in  a  FGM  layer  with 
E2/Ei=5  and  a2/ai=2  undergoing  a  uniform  temperature  change  T-To,  where  a  is  the 
thermal  expansion  coefficient.  Note  that  the  temperature  is  normalized  with  respect  to  Tq 
and  (T-To)/To=(T/To)-1.  Hence,  the  stress  <Tyy(x)  is  lienarly  dependent  on  (T/To)-l.  In  all 
cases  the  normalizing  stress  is  defined  by 

a„=i:,airo/(l-v)  .  (39) 

The  stress  ayy(x)  for  a  steady-state  heat  conduction  problem  in  a  FGM  layer  with 
E2/Ei=10,  a2/ai=2  and  X2/Xi=10  is  shown  in  Fig.  21,  where  A,(x)  is  the  coefficient  of  heat 
conduction.  In  these  examples  the  functions  E(x),  a(x)  and  X(x)  are  assumed  to  be  of  the 
form  given  by  (37).  The  normalized  stress  intensity  factors  for  an  edge  crack  in  a  FGM 
layer  (Fig.  13)  imder  thermally  induced  stresses  given  by  Figures  20  and  21  are  shown  in 
Figures  22  and  23,  respectively.  In  these  figures  the  normalizing  stress  intensity  factor  is 

CT  0  -v/ft  where  Oq  is  defined  by  (39). 

The  analytical  details  of  the  solution  of  the  surface  crack  problems  and  further 
results  may  be  found  in  [15*]  and  [16*]. 
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3.5  Stress  Relaxation  Due  to  Periodic  Surface  Cracking 

The  mud-flat  type  periodic  surface  cracking  of  ceramics  or  ceramic  coatings  under 
residual  or  thermal  stresses  seems  to  be  a  commonly  observed  mechanism  through  which 
the  stresses  or  energy  are  relieved.  Usually  such  materials  are  under  in-plane  biaxial 
strains  of  equal  magnitude.  An  important  mechanics  problem  here  is  the  determination  of 
the  length  parameter  describing  the  surface  crack  periodicity  in  terms  of  the  biaxial  strain 
So  and  the  fracture  touchness  A  highly  simplified  technique  used  in  practice  largely 
for  screening  purposes  is  an  energy  balance  approach  [30].  It  is  usually  assumed  that  in  a 
single  “cell”  formed  by  surface  cracks,  the  energy  available  for  fracture  is  provided  by  the 
strain  energy  within  the  cell,  whereas  the  energy  needed  to  create  the  fracture  surfaces  is 
the  energy  sink.  The  major  assumption  here  is  that,  as  a  result  of  cracking,  the  entire 
strain  energy  inside  the  volume  of  the  cell  is  released,  which,  in  piecewise  homogeneous 
materials,  leads  to  the  following  energy  balance  relation: 

,  \J^=AbW,  U^=Sb(bJ2)  ,  W  =  elE/(\-v),  (40) 

where  Ur  is  the  strain  energy  per  cell,  Uq  is  the  dissipated  energy  per  cell,  A  is  the 
surface  area  of  the  cell,  b  is  the  crack  depth,  W  is  the  strain  energy  density,  S  is  the 
peripheral  length  of  the  crack  for  a  single  cell,  and  bJ2  is  the  energy  required  to  create  a 
imit  fracture  surface.  If  the  cell  is  an  n-sided  regular  polygon  with  center  to  apex  distance 
c,  then 

A  =  —  sin(— )  ,  S  =  2nc  sin(— )  .  (41) 

2  n  n 

The  shape  of  the  polygon  or  n  must  be  such  that  for  a  unit  surface  area  UrAJq  is 
maximum.  For  the  purpose  of  determining  n  one  may  assume  that  A  is  “unity”.  Thus, 
eliminating  c,  for  given  Sq  and  bJ2  the  condition  UR/UD=max.  may  be  replaced  by 
A/S=max.,  from  which  it  would  follow  that 

(yl  /  5)  =  -(«  tan(— ))"*'^  =  max.  (42) 

2  n 

Technically,  (42)  gives  the  expected  result,  namely  n=oo  or  a  circular  cell.  However, 
since  the  plane  cannot  be  divided  continuously  into  circles  of  equal  radii,  physically  this 
result  is  not  acceptable.  The  admissible  regular  polygons  are  triangle,  square  and 
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hexagon.  Since  A/S  is  a  monotonously  increasing  function  of  n,  the  acceptable  solution 
of  (42)  would  be  n=6  or  a  regular  hexagon,  which  is  generally  the  observed  result  [30]. 

For  hexagonal  cells  we  have  A  =  3yf3c^  /  2 ,  S=6c  and  from  (40)  we  find 


CE 


2 

0 


2(1 -V) 
yf3E 


(43) 


where  c  is  the  length  of  the  side  of  the  hexagon.  Equation  (43)  gives  the  relationship 
between  the  cell  size  and  the  biaxial  strain  in  the  general  periodic  surface  cracking 
problem.  Note  that  in  this  simple  model  the  criterion  formd  is  independent  of  the  crack 
depth. 

Similar  results  may  be  obtained  if  we  assume  the  cracking  medium  (constrained 
single  layer  or  coating)  is  FGM.  For  hexagonal  surfaee  cracking  from  the  energy  balance 
equation  Ur=Ud  it  may  easily  be  shown  that 

b 

cEl=^i\bXx)dx)/i\Eix)dx/(\-v))  (44) 

v3  0  0 

where  and  E  are  the  fracture  toughness  and  Young’s  modulus  of  FGM  and  b  is  the 
crack  depth.  Note  that  in  FGMs  the  cell  size  vs.  biaxial  strain  relationship  foimd  is  now 
dependent  on  the  crack  depth.  It  is  possible  that  (44)  may  again  be  useful  for  screening 
purposes  if  one  assumes  the  crack  depth  to  be  equal  to  a  certain  percentage  of  the  FGM 
coating  thickness  h. 

The  preceding  technique  is  very  highly  approximate.  Since  the  corresponding 
three-dimensional  elasticity  problem  for  the  periodic  surface  cracks  is  analytically 
intractable,  to  shed  some  light  on  the  subject  one  may  consider  a  simple  two-dimensional 
problem.  The  simplest  such  problem  of  periodic  surface  cracking  of  a  FGM  coating 
bonded  to  a  homogeneous  substrate  vmder  antiplane  shear  loading  was  recently 
considered  in  [17*].  In  this  study  it  was  assumed  that  the  substrate  is  a  half  space  having 
the  shear  modulus  \X2  and  the  shear  modulus  of  the  FGM  is  given  by 


^i(^)  =  lii(0)exp(Px)  ,  p=ilog(-^) 

h  p,(0) 


(45) 
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where  )ai(0)  is  the  shear  modulus  of  FGM  on  the  surface,  x=0,  and  h  is  the  thickness  of 
the  coating.  It  is  further  assumed  that  the  coating  contains  periodic  surface  cracks  with 
the  period  2c  and  depth  b,  and  the  medium  is  subjected  to  a  constant  strain 
2e  ^  (x,+oo)  =  y  0  through  fixed  grips  (Fig.  24). 

The  results  of  an  example  showing  the  stress  relaxation  on  the  surface  of  the 
coating  as  a  result  of  periodic  cracking  are  shown  in  Figures  25  and  26.  Note  that  as 
c/h->^0  the  surface  stress  criyz(0,y)  would  tend  to  zero  everywhere  and,  as  c/h-)-oo, 
aiyz(0,c)  would  approach  Tq,  the  surface  shear  stress  for  the  uncracked  coating.  We  also 
note  that  the  previous  simple  periodic  cracking  model  assumes  the  stress  to  be  zero 
within  the  entire  cell,  that  is,  aiyz(0,y)=0,  0<y<c.  Figures  25  and  26  show  that  this  is 
hardly  the  case. 

Another  interesting  result  of  this  study  is  shown  in  Figures  27  and  28.  Here  Vj  is 
the  actual  strain  energy  released  per  unit  area  of  the  coating  as  a  result  of  periodic 
cracking.  Vj  is  calculated  from  the  elasticity  solution  for  the  cracked  medium.  Ui  is  the 
total  strain  energy  contained  in  the  volume  1x1  xb  and  is  calculated  from 

6 

U,=UJlc  ,  U,=cyl  Jp(x)n^r  .  (46) 

0 

Clearly,  the  physics  of  the  problem  requires  that 

limr,  =C/,  .  (47) 

c->0 

Figures  27  and  28  show  that  the  assumption  made  for  the  previous  simple  model  to  the 
effect  that  the  entire  strain  energy  contained  within  the  cell  is  relieved  and  is  transformed 
into  fracture  energy  may  be  justified  only  for  very  small  values  of  c/h  and,  in  general, 
only  a  certain  percentage  of  the  total  energy  is  available  to  overcome  the  fracture  energy. 

Figures  29  and  30  show  the  way  one  may  use  the  basic  results  giving  the  stress 
intensity  factor  as  a  function  of  the  crack  depth  in  fracture  stability  studies  in  FGM 
coatings.  Assuming  that  near  the  surface  the  FGM  is  100%  ceramic  and  near  the 
interface  100%  metal,  in  this  example  the  mode  III  fracture  resistance  parameter  k3(, 
obtained  from 

=  (48) 

21 


is  taken  to  be  an  inereasing  function  of  x  as  shown  in  the  figures.  Consequently,  the 
FGM  coating  would  exhibit  an  R-curve  (crack-extension-resistance-curve)  behavior. 
That  is,  since  the  fracture  toughness  increases  with  increasing  crack  length,  the  load 
amplitude  must  be  increased  continuously  to  keep  the  crack  growing.  Thus,  the  threshold 
value  =  Pi(0)y  o^.  and  the  critical  value  t  ^  =  lii(0)y  may  be  obtained  by  requiring 
that  the  corresponding  k3(b)  curves  envelope  k3c(x)  from  below  and  from  above, 
respectively.  For  T7-<To<T^the  crack  would  grow  in  an  imstable  manner  in  the  region 
for  which  k3>k3c  and  becomes  arrested  when  k3^3(.. 

Extensive  results  on  the  periodic  cracking  of  FGM  coatings  under  Mode  III 
loading  conditions  may  be  found  in  [17*].  The  plane  strain  problem  of  periodic  cracking 
for  a  homogeneous  coating  and  for  a  single  homogeneous  layer  subjected  to  membrane 
loading  or  bending  was  considered  in  [18*]. 

3.6  End  Effects  on  Stress  Concentration  and  Edge  Debonding 

Generally  the  stress-free  ends  in  bonded  materials  are  locations  of  high  stress 
concentrations  and  potential  debonding  fracture.  In  bonded  dissimilar  homogeneous 
materials  the  point  at  which  the  interface  intersects  the  free  boimdary  (or  the  apex  of  two 
90  degree  bonded  wedges)  is,  in  fact,  a  point  of  singularity  near  which  the  stress  state  is 
given  by  [31] 

=  -  ('.>  =  *..>’)  .  0<p<l/2  ,  (49) 

when  (r,0)  are  the  polar  coordinates,  p  and  Fy  depend  on  the  bimaterial  constants  and  K  is 
a  measure  of  the  load  amplitude  or  stress  intensity.  For  p  to  be  positive  the  material 
properties  need  to  be  discontinuous  across  the  interface.  In  FGM  coatings,  since  the 
material  properties  are  made  continuous  through  composition  grading,  it  can  be  shown 
that  the  singularity  P  becomes  zero  and,  consequently,  the  stresses  become  finite  [32]. 

Figure  31  describes  the  geometry  of  a  specimen  considered  as  an  example.  The 
composite  medium  is  assumed  to  undergo  a  homogeneous  temperature  change  AT.  The 
problem  is  one  of  plane  strain.  The  substrate  is  a  nickel-based  superalloy  (Rene  41)  and 
the  coating  is  either  a  piecewise  homogeneous  or  a  functionally  graded  metal/ceramic 
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layer,  the  ceramic  component  being  partially  stabilized  zirconia  (PSZ).  The 
thermomechanical  constants  of  the  two  materials  are 

Es-219.7  GPa ,  V3=0.3  ,  as-1.67xlO‘^°K'‘ , 

Ec=151  GPa ,  v,=0.3  ,  ae=10‘^°K 

where  the  subscripts  s  and  c  refer  to  the  substrate  and  the  ceramic,  respectively.  The 
stepwise  variation  of  the  material  properties  in  piecewise  homogeneous  coatings  is  shown 
in  Fig.  32.  Some  calculated  results  are  shown  in  Figures  33-37.  For  the  FGM  coating 
used  in  the  example,  the  modulus  variation  is  given  by 

Ej, ,  0<y<0.0125  m. 

E(y)  =  { 

00145-  V 

Ee+(E3-Ec)(  ,  0.0125<y<0.0145  m. 

The  same  expression  is  ued  for  the  thermal  expansion  coefficient,  a.  Figure  38  shows  the 
thickness  variation  of  the  stiffness  E(y)  (or  the  thermal  expansion  coefficient  a(y)  for  the 
three  FGMs  considered,  namely  the  metal-rich,  linear  and  ceramic  rich  compositions. 

Figures  33  and  34  show  the  interface  stresses  CTyy(x,h2)  and  axy(x,h2)  for  the 
piecewise  homogeneous  coating.  Note  that  since  the  material  properties  are 
discontinuous  for  y=h2,  the  interface  stresses  become  imboimded  for  y=h2,  x->L  The 
corresponding  results  for  the  FGM  coatings  are  shown  in  Figures  35  and  36.  It  may  be 
observed  that  as  a  result  of  material  property  smoothing  (or  eliminating  the  property 
discontinuities),  the  stress  singularities  are  eliminated  and  at  the  point  y=h2,  x=f  CTyy(x,h2) 
becomes  bounded  and  CTxy(x,h2)  becomes  zero.  One  may  also  observe  that,  except  for 
values  of  x  near  the  end  (x={),  the  behavior  of  interface  stresses  for  piecewise 
homogeneous  and  FGM  coatings  are  quite  similar.  Figure  37  shows  the  variation  of  the 
tensile  stress  Cyy((,y)  at  the  ends  of  the  specimen  for  0<y<hi+h2.  Note  that  at  the  interface 
y=h2  there  is  a  stress  concentration  and  the  stress  concentration  factor  Cyy/ao  decreases 
with  increasing  metal  content  of  the  coating.  Further  results  on  this  problem  may  be 
found  in  [8*]. 

The  effect  of  material  orthotropy  on  the  crack  driving  force  for  homogeneous 
debonding  layers  was  considered  in  [19*]  and  [20*].  In  [19*]  the  plane  strain  problem  of 
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debonding  crack  initiating  from  a  free  end  was  studied.  The  problem  of  a  T-shaped  crack 
simulating  the  debonding  process  that  is  initiated  from  a  surface  crack  terminating  at  the 
coating/substrate  interface  was  considered  in  [20  ]. 

3.7  Some  Concluding  Remarks 

From  the  viewpoint  of  fracture  mechanics  the  “functionally  graded  materials” 
seem  to  offer  certain  advantages  among  which  one  may  mention  the  following: 

(i)  By  eliminating  the  discontinuity  in  material  property  distributions,  the 
mathematical  anomalies  regarding  the  crack  tip  stress  oscillations  for  the  interface  cracks 
and  the  non-square-root  singularities  for  the  cracks  intersecting  the  interfaces  are  also 
eliminated.  In  practice  the  importance  of  this  result  lies  in  the  fact  that  in  FGMs  one  can 
now  use  the  crack  tip  finite  element  modeling  developed  for  the  ordinary  square-root 
singularity  and  apply  the  methods  of  the  energy  balance-based  theories  of  the 
conventional  fracture  mechanics. 

(ii)  Application  of  FGMs  as  interfacial  zones  in  joining  generally 
incompatible  materials  would  greatly  improve  the  bonding  strength. 

(iii)  Use  of  FGMs  as  coatings  and  interfacial  zones  would  reduce  the 
magnitude  of  the  residual  and  thermal  stresses. 

(iv)  Use  of  FGM  coatings  and  interfaces  would  eliminate  the  stress 
singularities  at  the  points  of  intersection  of  interfaces  and  stress-free  ends  in  bonded 
materials. 

(v)  Replacing  homogeneous  coatings  by  FGM  layers  would  both  enhance  the 
bonding  strength  and  reduce  the  crack  driving  force. 
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Figure  1:  Notation  for  a  plane  crack  in  a  nonhomogeneous  medium. 


ki2.  kli  ~5  (t-x)  kts,  kli  ~  >09  M  ><i2.  ><21  ~  Cont. 


Figure  3:  Singular  behavior  of  the  irregular  kernels  for  an  interface  crack  in  bonded 
dissimilar  materials. 


Fig.  4:  Geometiy  of  the  interface  crack  for  a  functionally  graded  coating  bonded  to 
homogeneous  substrate. 


ki/ko  0.8  1 


-0.2 


-3 


1  2 


—I 

3 


Fig.  9  The  influence  of  *2  /  a  on  fcj ,  loading:  x  „ ,  v=0.3,  /  a  =  100,  (1)  ^  ° 

(2)  /i2  /  ^  ♦  (^)  ^2  !  0.25 . 


Fig.  11  The  effect  of  the  coating  ttuckness  hi/  a  and  the  nonhomogeneity  constant  y  on 
the  probable  crack  growth  direction  G,,,  loading;  Cg,  v=0.3,  hi/a  =  100,  (1) 
/i2 /a  =  10  100,  (2)  hi!  0  =  2,0)  h2la  =  \A^)  *2  /  a  =  0.5 .  (5)  h2/a  =  025. 

(The  angle  6  ^  is  measured  at  the  crack  tip  a  from  the  positive  x  axis.) 


Fig.  12  Normalized  stress  -intensity  factors  for  a  plane  crack 
terminating  at  the  interface  between  a  homogeneous  medium  and 
a  FGM  half  space. 
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Fig.  14  The  stress  distribution  in  the  uncracked  layer  under  fix 

CTfl  =  E2eg/(1  -v^). 
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Fig.  17  Stress  intensity  factor  for  an  edge  cra< 
loading,  a  =  0,  (Tq  =  EjeQ(l  — 
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Fig.  18  Stress  intensity  factor  for  an  edge  era' 

loading,  a  =  0,  =  N/h. 
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Fig.  20  Thermal  stress  distribution  in  a  FGM  layer  with 
E2/Ei=5,  02/011=2,  undergoing  a  uniform  temperature  change, 
Ti=T2=T,  T<Tq. 
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Fig.  21  Thermal  stress  distribution  in  a  FGM  layer  with 
E2/E2=10,  a2/a2=2,  X2/^l=10,  under  steady-state  heat 
conduction;  T(h)=T2=0.5TQ,  T(0)=Ti  variable. 
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Fig,  23  Normalized  stress  intensity  factor  for  a  surface  crack  in 
a  FGM  layer  under  the  stress  state  shown  in  Fig.  21  (E2/E2^=10, 

X2/Xi=10)  . 
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Fig.  25  The  stress  <7iy,(0,j/)  on  the  surface  of  a  periodically  cracked  FGM 
coating  bonded  to  a  homogeneous  medium  under  remote  constant  shear 
strain  7o;  =  0.5,  b/h  =  0.1,  Tq  =  7o/ii(0). 
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Fig.  26  Same  as  Figure25  =  0.5,  b/h  =  0.5,  Tq  =  7o^(0). 
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Fig.  27  Released  strain  energy  per  unit  surface  area  of  a  periodically 
cracked  FGM  coating  bonded  to  a  homogeneous  substrate  under  constant 
strain  loading  70;  =  1.0,  Ui  =  (l/2))Ui(0)7o(e^^  — 

1 .0 

0.8 

0.6 

0.4 

0.2 

0.0 

0.0  0.5  1.0  1.5  2.0  2.5  3.0 

c/h 

Fig.  28  Same  as  Figure  27  ^A=s  — 1.0,  — l)//S. 


Variation  of  Stiffness,  E  (GPa) 
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Fig.  31  The  geometry  of  the  coating/substrate  considered 
as  an  example. 
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Fig.  32  Stepwise  distribution  of  the  elastic  modulus 
for  a  single,  two  and  four  layer  homogeneous  coatings. 
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Figure  33  Variation  of  normalized  along  the  1st  interface  of  laminated  homogeneous 
coating;  Gg  =  E^a^AT. 


Figure  34  Variation  of  normalized  along  the  1st  interface  of  laminated  homogeneous 
coatings;  Cg  =  E^a^AT. 
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Figure  37  Variation  of  normalized  <Sy  along  the  end  x  =  i  —  0.Q2m  for  FGM  coatings; 
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Fig*  38  Distribution  of  the  elastic  modulus  E(y)  for  the  FGM  coating. 


